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THE RED-SHIFT— MAGNITUDE RELATION AND OBSERVATIONAL DATA 


By Windsor L. Sherman 
Langley Research Center 

SUMMARY 

The red-shift — magnitude relation for uniform, zero-pressure models of the uni- 
verse is one of the stronger relationships connecting theory and observation in cosmology. 
This relationship has been used to determine models for the universe from the observed 
apparent magnitudes and red shifts of galaxies. The problems associated with the use of 
the red-shift — magnitude relation to determine a model universe are discussed. In par- 
ticular, the effects of evolution and selection on the use of the red-shift — magnitude rela- 
tion are considered. Selection is the tendency of observers to select intrinsically brighter 
sources for observation as the distance to the sources increases. The results indicate 
that the red-shift — magnitude relation should be used with caution when determining a 
model for the universe because the variation of the intrinsic luminosity of sources must 
be taken into account. 


INTRODUCTION 

The theory of uniform models of the universe based on general relativity provides 
several connections between theory and observation for models in which the pressure is 
equal to zero. Sherman (ref. 1), McVittie (ref. 2), Solheim (ref. 3), and Sherman and 
Wallace (ref. 4) have studied the connections between theory and observation and derived 
exact and simplified relationships for the analysis of observational data. Although the 
use of this connection between theory and observation has been illustrated, there has been 
little discussion of the problems of calculation and interpretation of the results. 

In this report some of the problems connected with the use of the red-shift — 
magnitude relation to obtain a model of the universe from observational data are dis- 
cussed. Consideration is given to the effects of evolution and selection on the model. 
Selection is the tendency of observers to select intrinsically brighter sources for obser- 
vation as the distance to the sources increases. The results obtained from the red-shift — 
magnitude relation are compared and interpreted with respect to models based on other 
observational data and different methods of model determination. The symbols used in 
the analysis are defined in appendix A. 


THE RED-SHIFT— MAGNITUDE RELATION 


The most general form of the red-shift — magnitude relation for zero-pressure 
relativistic models of the universe is (ref. 4) 


where 


m - K = 5 log 10 [Ro(l + z)S(cl>)J + M + AM - 5 
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The variable co is given by an elliptic integral of the first kind as 


o> = \J\3a 0 - q D - l| 2a 0 z 3 + (3a Q + q 0 + l)z 2 + 2(q 0 + l)z + 1 


- 1/2 


dz 
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The correction for evolutionary effects AM is given by 

AM = -M(t 0 - t) 

and the time of light travel t G - t is 


t 0 ~ t = H 0 
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In the foregoing equations, m is the apparent magnitude, K is the red-shift correction, 
R 0 is the present value of the scale factor, z is the red shift, M is the absolute mag- 
nitude, c is the speed of light in a vacuum, H 0 is the Hubble parameter, o 0 is the 
density parameter, q 0 is the deceleration parameter, M is the rate of change of abso- 
lute magnitude, t 0 is the present time, and t is any other time. When the red-shift — 
magnitude relation is used to obtain a model of the universe, the unknown parameters are 
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H 0 , (Tq, q 0 , and M. However, use of the method of differential corrections showed that 
M and H 0 are too interrelated and highly correlated to permit these two parameters 
to be determined. The parameters H 0 and M are usually combined into a single 
parameter, and when this is done, equation (1) is written as 


m - K = 5 log 1Q At , 2 . S(w) + AM + C (7) 

\j |3°b - q Q - 1| 

where 

C = 5 log 10 jA + M - 5 (8) 

Combining H 0 and M into a single parameter means that only C, o 0 , and q Q need 
be determined from the observed data. After C, cr 0 , and q Q have been obtained from 
observational data, they are used to specify a model for the universe in terms of the 
curvature constant k and the cosmical constant A. The constant k is given by equa- 
tions (4) and 

A = 3H 0 ( °o “ %) 

The parameter C (eq. (8)) is made up of the Hubble parameter H 0 , the speed of 
light in a vacuum c, and the absolute magnitude of the source M. The Hubble param- 
eter which is determined from the observation of nearby galaxies is assumed to be a con- 
stant. The speed of light in a vacuum is also a constant. Thus, the absolute magnitude 
is the only quantity that can vary in the expression for C. The absolute magnitude is a 
function of the physical process that is supplying energy for the source, the state of the 
process, and the size and composition of the source. Because none of these factors are 
a function of the position of the source, the absolute magnitude is independent of the model 
universe. In a complicated structure like a galaxy where the emission is the integrated 
output of many stars and nebulae, the absolute magnitude depends on the number of emis- 
sion sources and the state of the energy process of each source. Inasmuch as the size of 
a galaxy is a rough indication of the number of sources, the absolute magnitude of gal- 
axies varies. Thus, if n sources are being used to determine a model for the universe, 
it may be necessary to determine n values for C in addition to the parameters cr 0 
and q Q . Because it is not possible to determine n + 2 values of the parameters C, 
a Q , and q Q , an assumption must be introduced to render the problem determinant. Inas- 
much as ct q and q G are the parameters that determine the model, the usual assump- 
tion introduced is that all sources in the set of data being used have the same value of C. 
This assumption is equivalent to saying that all sources have the same absolute magnitude. 


ANALYSIS OF OBSERVATIONAL DATA BY USE OF THE 


RED-SHIFT— MAGNITUDE RELATION 


The form of the red-shift — magnitude relation given in equation (7) and the method 
of differential corrections were used to obtain the parameters a Q , q 0 , and C from the 
observational data for clusters of galaxies. The best values of the parameters cr 0 , q Q , 

and C were obtained when the sum of the squares S (m n - m ca i)^ was a minimum. 

i=l 1 

The values of ct 0 , q Q , C, and the root mean square (rms) on the magnitude residuals 

\ EK - m cal)f 

were obtained. The rms is given by y where m c is the corrected 

observed apparent magnitude, m ca i * s calculated apparent magnitude obtained from 
the red-shift — magnitude relation, and n is the number of sources in the group being 
analyzed. The observational data used to determine a Q , q Q , and C are presented in 
table I, discussed in appendix B, and plotted in figure 1. 

TABLE I.- OBSERVATIONAL DATA FOR CLUSTERS OF GALAXIES 


Cluster 

Red shift, 
z 

Corrected apparent 
magnitude , 
m c 

Virgo 

0.004 

9.27 

0316 + 4121 

.018 

12.51 

1257 + 2812 

.022 

12.91 

1603 + 1755 

.036 

14.12 

2308 + 0720 

.043 

14.78 

2322 + 1425 

.044 

15.04 

1145 + 5559* 

.052 

15.71 

0106 - 1536 

.053 

15.21 

1024 + 1039 

.065 

15.88 

1239 + 1852* 

.072 

15.22 

1520 + 2754* 

.072 

15.86 

0705 + 3506 

.078 

16.26 

1431 + 3146* 

.131 

17.31 

1055 + 5702 

.134 

17.20 

0025 + 2223 

.159 

17.39 

0138 + 1840 

.173 

17.16 

0925 + 2044 

.192 

17.56 

0855 + 0321 

.202 

17.84 

0024 + 1654 

.290 

18.88 

1448 + 2617* 

.350 

18.68 

1410 + 5224 

.460 

19.48 


*Does not fall nicely on a simple curve. 
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Figure 2 shows the results 
obtained when the data were pro- 
cessed to determine ct 0 , q 0 , and 
C. In this calculation, M was set 
equal to zero and it was assumed 
that all sources in the group had the 
same value of the parameter C. 

The values of cr 0 and q Q found 
for this set of data give a curvature 
constant k of +1 and a cosmical 
constant A that is greater than 
zero. In addition to the o Q , q Q , 
and C values obtained from pro- 
cessing the data, the present value 
of the scale factor R 0 , the time 
since the beginning of expansion T, 
the variation of R/Ro with T, the 
variation of t Q - t with z, the 
time of light travel, and the varia- 
tion of m with log^Q z are a l so 
given in figure 2. The equations 
and method for determining T, 

R 0 , and R/R 0 are given in refer- 
ence 4. The information presented 
in figure 2 gives a complete 
description of the model universe. 
The values of a Q and q Q of this 
model agree well with the values of 
a Q and q Q of the Solheim model 
(ref. 3). The values of a 0 and q 0 
found by Solheim were 4.53 and 
-0.06, respectively. The small dif- 
ferences between cr 0 and q Q 
given in figure 2 and those obtained 
by Solheim probably arise from the 
omission of the Virgo cluster from 
Solheim' s data. The model given 
in figure 2 was obtained from 
observational data and should 



Figure 1.- Observational data for clusters of galaxies 
(from table I). 


Computed parameters 




Figure 2 - Model for the universe obtained when the observational data 
from table I were analyzed by using the red-shift— magnitude relation 
with parameter C assumed to be a constant and M = 0. 
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represent the observed universe if general relativity is the theory of gravitation that 
applies to the universe. It is possible to obtain independent estimates of the density 
parameter and the time since the beginning of expansion through the study of the Galaxy 
and clusters of galaxies. Because the data from which the model obtained by means of 
the red-shift — magnitude relation come from the same universe as the independently 
determined values of a Q and T , reasonable agreement should exist between the two 
sets of values. Lastly, the model should have a feature or features that can be used to 
account for the 3° K blackbody radiation that appears as an isotropically distributed back- 
ground in the universe. 

Comparison of Red-Shift — Magnitude Model Parameters With 
Independently Determined Values 

The independently determined values of density parameter and time since the 
beginning of expansion are based on estimates of the present density of the universe and 
the age of the Galaxy. These data and the 3° K blackbody radiation are discussed in 
appendix B. 

Density .- For the model given in figure 2, 2.503 S a Q £ 6.277 with the most likely 
value of the density parameter being 4.39. The independent determination gives 
0.0045 s cr 0 s 0.45 with the most likely value being 0.045. The density parameter inter- 
val for the model of figure 2 is completely different from the density parameter interval 
based on observational results, and the most likely values differ by a factor of 97.6. 

Time since beginning of expansion .- In a model, like that shown in figure 2, where 

the expansion started at zero radius, the initial state was extremely compact and it is 

generally hypothesized that the galaxies formed very soon after expansion started. This 

hypothesis means that age of all galaxies is about the same; thus, the time since the 

beginning of expansion should be greater than or at least equal to the age of the Galaxy. 

The time since the beginning of expansion, based on the variance of a Q and q 0 , falls 

in the interval 3.65 X 10^ st ^ 6.86 x 10^ years with the most likely value being 

4.48 x 10^ years. These times are for H 0 = 3.24 x 10"-*-® sec"-*-. The age of the Galaxy 

A a falls in the interval 1.2 x 10*-® ^ Ap- = 2.0 X 10-*-® years. The most likely value of 
° in 5 

the age of the Galaxy is 1.5 x 10 x ^ years. Comparison of these data show that the time 
interval for the age of the Galaxy and the time interval for the time since the beginning 
of expansion are completely different. The most likely age of the Galaxy is about 3^ 
times the most likely time since the beginning of expansion. In a model universe such 
as that shown in figure 2, a galactic age greater than T is not admissible. 

The 3° K blackbody radiation .- An isotropically distributed 3° K blackbody radia- 
tion has been observed in the universe. The present hypothesis concerning this radiation 
is that it is the remnant of the initial state of the universe. In this initial state the 
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universe was an extremely compact, hot, dense body that expanded into the present 
observed universe, the so-called big bang model. The model given in figure 2 has a zero 
radius (radius is proportional to R) at zero time. This feature would give an extremely 
compact initial state and, thus, the conditions for the development of the 3° K blackbody 
radiation are found in this model. 

D iscussion of t he comparison .- Except for the 3° K blackbody radiation, the model 
universe presented in figure 2 does not have a density parameter and time scale that are 
compatible with independently determined values of these quantities. This lack of com- 
patibility could be due to (1) use of an incorrect theory of gravitation, (2) too much simpli- 
fication of the basic model, (3) incorrect data, and (4) incorrect use of the data leaving 
out the effects of evolution and selection or use of an incorrect assumption when the data 
were analyzed. Because there is insufficient evidence for either item (1) or item (2) and 
saying that the data are incorrect as in item (3) amounts to questioning the validity of a 
large amount of astronomical data without a good basis, it was decided to follow the course 
set up in item (4). 


Uncertainties in Red-Shift — Magnitude Relation 

Besides the assumptions made in the derivation of the red-shift — magnitude rela- 
tion, the most important assumption is that the parameter C is a constant. This 
assumption implies that all galaxies have the same absolute magnitude. The absolute 
magnitude of all galaxies is not the same. For instance, the galaxies marked with an 
asterick in table I could not have the same absolute magnitude as the adjacent galaxies. 
See reference 5 (p. 272) for absolute magnitudes of galaxies determined from Cepheid 
variables. Because all galaxies do not have the same absolute magnitude, a study was 
made to determine the effect of assuming constant absolute magnitude on the model when 
the actual absolute magnitudes were not the same. The study was made with the red- 
shift — magnitude relation and a synthetic set of data. In the synthetic data each source 
had a different value of C and this parameter was assumed to be a constant in the red- 
shift — magnitude relation. It was found that when intrinsic luminosity increased with 
distance, the assumption of constant intrinsic luminosity in the red-shift — magnitude 
relation caused the extracted value of a Q to be much greater than the value of cr 0 used 
in the construction of the synthetic data; that is, a high-density universe was predicted 
instead of the low-density universe used to obtain the data. The evolution of galaxies and 
selection, both observational and instrumental, could produce a set of data in which the 
distant sources were brighter than the near sources. 

Effect of evolution on model for the universe .- The evolutionary effect is based on 
the change of the luminosity of galaxies with time and the time of light travel from the 
source to the observer. Evolution has been discussed by several authors (for instance, 



refs. 2 and 6). A value of M of 0.3 per 10 s * years was selected to be a value of M 
large enough to illustrate the effect of evolution and small enough not to be unrealistic. 
The data were analyzed to determine a 0 and q Q with M equal to the foregoing value. 
The results of this calculation are presented in figure 3. Comparison of the data pre- 
sented in figure 3 with those presented in figure 2 shows that evolution reduces both a 0 
and q 0 . The density parameter now lies in the interval 0.22 = cr 0 = 3.0 with the most 
likely value being 1.61. The change in ct 0 reduces the discrepancy with the indepen- 
dently obtained value. The time since the beginning of expansion now lies in the interval 
4.79 x 10^ years with the most likely value being 6.59 x 10^ years, an increase in T 
that reduces the discrepancy with the age of the galaxy. These times are for 
H 0 = 3.24 x 10 - -*-® sec - l. The model of figure 3 has a zero radius at t = 0; thus, the 
3° K blackbody radiation can be accounted for in this model. The introduction of evolu- 
tion caused changes in ct 0 and q 0 which apparently moved the model toward the 
observed universe; this indicates the importance of accounting for evolution when deter- 
mining a model for the universe. This result agrees with the result of reference 6 con- 
cerning the importance of evolution. Even though the changes in ct 0 and q Q improved 
the model, this change is not reflected in the rms of the magnitude residuals. The rms 


Computed parameters 



Figure 3.- Model for the universe obtained when the observational data 
from table I were analyzed by using the red-shift— magnitude relation 
with M = 0.3 per 10 9 years. 
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is the same as that for the model given in figure 2. The lack of change in the rms is 
more than likely caused by a high dependence of the evolutionary effect on ct 0 and q D . 
This dependence of the evolutionary effect on a Q and q 0 apparently absorbs the effects 
of evolution on the rms. 

Selection effect .- Selection is a complicated process and involves such factors as 
the instruments, observers, and available sources. Both the size and physical state of 
the available sources enter into the problem of selection. In observational cosmology, 
such factors as instrument sensitivity, recording time, reading of survey plates, and 
available sources influence the observer to select the brightest of the available sources 
for study. At large distances only the brighter sources are available; therefore, there 
is considerable chance that the distance sources will be intrinsically brighter than the 
nearer sources. A discussion of selection and how it can affect the model for the uni- 
verse is presented in reference 7; however, the effects of selection have not yet been 
detected in the observational data. 

A study of observational data was made to determine whether any trends in the data 
might indicate the presence of a selection effect. The apparent magnitudes and red shifts 
of 223 galaxies, principally from reference 8 were separated into groups as shown in 
table II. In order to reduce and if possible to eliminate random errors, the data for the 
galaxies in each group were averaged by determining the mean apparent magnitude and 
mean red shift for the group. The mean apparent magnitude is given by 

Zw 

i=l 



and the mean red shift is given by 


lo §10 z k 

where n indicates the number of galaxies in the group and k refers to the group. In 
addition, the variance was also determined, that is, 

n 

Z( m l-“k ) 2 
a 2 = t± 

n 

The results of the calculations are also given in table II. 


n 

2 1Og 10 Z i 

i=l 

n 
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TABLE II.- MEAN APPARENT MAGNITUDE, VARIANCE, AND 
RED SHIFT OF 223 GALAXIES 


Red-shift interval 

m 

o2 

l°gl0 cz 

No. of 
sources 

0.0133 = z < 0.0167 

13.64 

0.81 

3.65350 

62 

0.0167 s z < 0.025 

14.28 

.75 

3.77759 

64 

0.025 Sz< 0.0333 

14.49 

.70 

3.93171 

18 

0.0333 Izs 0.05 

15.51 

.32 

4.07583 

16 

0.05 s z so. 0667 

15.93 

.70 

4.22897 

11 

0.0667 s z s 0.1 

16.96 

1.12 

4.38533 

18 

0.1 <z so. 1833 

18.28 

CD 

t> 

4.63838 

19 

0.1833 s z 

19.11 

.41 

4.83982 

15 


The mean magnitudes given in table II were corrected for the effects of red shift 
(ref. 8). The corrected mean magnitude data are plotted in figure 4 with logjp z as 
the ordinate. In addition to mj Cc the variation of m with logjQ z for a model uni- 
verse in which o 0 = q Q = 1.0, the linear model, is also shown in the figure. The bars 
through the symbols for indicate the la variation for the computed mean at these 

points. For small values of z, the curve of m plotted against logjg z for the linear 
model lies within the error bars of m^. However, as log^o z (and, hence, distance) 
increases, the separation of mj C( . caused by the fact that mj^ is brighter than m at 

the same logjQ z > keeps the model line out 
of the error bars. The difference between 
mj (c and the model line is made up of differ- 
ences due to the model (that is, a 0 and q Q 
different from 1.0), evolution, and selection. 
The dashed curve in figure 4 shows the effect 
of model and evolution in the interval 
-1.8236 s logjQ z = -0-6373. The upper limit 
corresponds to a red shift of about 0.23. 

This curve is for the model given in figure 3 
and is considered typical for a model with 
evolution in the interval for log^g z - There 
is not an inconsiderable difference between 
the dashed curve and the mk c points, and 
as the upper limit is approached, the error 
bars no longer touch the model curve. These 
differences can be attributed to selection. In 
fact, figure 4 shows that if selection is 
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Figure 4.- Comparison of corrected mean apparent 
magnitude with la error bars and models for the 
universe. 
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present, it may be greater than or equal to the model increment; this means that with 
present data the red-shift — magnitude relation should not be used to determine a model 
for the universe unless the parameter C can be determined by independent means for 
each source. 


An Alternate Approach to Determination of a Q and q D 

Several methods have been studied for making an independent estimate of the param- 
eter C. The most satisfactory method depends on the noise in the recorded data and 
estimates of the possible range of values of cr 0 and q Q . This method, discussed in 
appendix C, was used in conjunction with data given in table I to estimate parameter C. 
The results of the estimation of C are given in table III. The method used to estimate 
the parameter C gives rise to many possible combinations. At least 272 combinations 
involving C are possible but only one is given in table HI. The combination presented 
in table III was selected so that the most distant sources were intrinsically brighter than 
the nearer sources. To make use of these values of C, m-Q, the apparent magnitude 
for the model a Q = q Q = 1.0, was determined for each source by use of 

m n = 5 log 10 z + C (9) 

and was then subtracted from the observed apparent magnitude. This procedure gives 
the incremental apparent magnitude Am. 

TABLE HI.- POSSIBLE VALUES OF THE PARAMETER C AND OF Am FOR 
OBSERVATIONAL DATA FOR CLUSTERS OF GALAXIES FROM TABLE I 


Cluster 

z 

5 log 10 z 

m c 

C 

Am 

Virgo 

0.004 

-11.990 

9.27 

21.26 

0 

0316 + 4121 

.018 

-8.725 

12.51 

21.26 

-.025 

1257 + 2812 

.022 

-8.290 

12.91 

21.26 

-.06 

1603 + 1755 

.036 

-7.220 

14.12 

21.26 

.08 

2308 + 0720 

.043 

-6.835 

14.78 

21.66 

-.045 

2322 + 1425 

.044 

-6.785 

15.04 

21.66 

.165 

1145 + 5559 

.052 

-6.420 

15.71 

22.06 

.07 

0106 - 1536 

.053 

-6.380 

15.21 

21.66 

-.07 

1024 + 1039 

.065 

-5.934 

15.88 

21.66 

.154 

1239 + 1852 

.072 

-5.715 

15.22 

20.86 

.075 

1520 + 2754 

.072 

-5.715 

15.86 

21.66 

-.085 

0705 + 3506 

.078 

-5.540 

16.26 

21.66 

.14 

1431 + 3146 

.131 

-4.415 

17.31 

21.66 

.065 

1055 + 5702 

.134 

-4.365 

17.20 

21.66 

-.095 

0025 + 2223 

.159 

-3.995 

17.39 

21.26 

.125 

0138 + 1840 

.173 

-3.810 

17.16 

20.86 

.11 

0925 + 2044 

.192 

-3.585 

17.56 

20.86 

.285 

0855 + 0321 

.202 

-3.475 

17.84 

21.26 

.055 

0024 + 1654 

.290 

-2.690 

18.88 

21.26 

.31 

1448 + 2617 

.350 

-2.280 

18.68 

20.86 

.10 

1410 + 5224 

.460 

-1.685 

19.48 

20.86 

.305 



The incremental apparent magnitude Am and the red shift for each source were 
then used in the equation 

Am = 5 log 10 - — A + Z == S(u>) (10) 

z^|3a 0 - q Q - 1| 

to determine a 0 and q Q . This equation was developed in the paper for which the 
abstract is presented in reference 9. The results of this calculation are shown in fig- 
ure 5. It was found that ct 0 falls in the interval 0 ^ cr 0 g 0.779 with a most likely 
value of 0.217. The most likely value of a Q based on observation now falls in this 
interval. The results in figure 5 show that q Q falls in the interval 0.0212 = q 0 = 0.782 
with the most likely value being 0.402. Compared with the values for the models in fig- 
ures 2 and 3, a Q has decreased and q 0 has increased. The decrease in cr 0 is in the 
correct direction to increase the time since the beginning of expansion. This parameter 
now lies in the interval 5.92 x 10^ $T i 9.71 x 10^ years with the most likely value of 
7.29 x 10^ years, which moves the time since the beginning of expansion closer to the age 
of the Galaxy. These times are for H 0 = 3.24 x 10-18 

sec - -*-. In the model presented in 

figure 5, the rms of the magnitude residuals has improved by a factor of more than 3 


Computed parameters 


m 91 7+ 



Figure 5.- Model for the universe obtained when the data in table III were 
analyzed by using the incremental red-shift— magnitude relation (eq. (10)). 
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over the previous models. The fact that this model has a zero radius at t = 0 means 
that it is simple to account for the 3° K blackbody radiation in this model. These results 
indicate that the inconsistency between the independently determined values of the density 
parameter and the age of the universe is much reduced when the parameter C is inde- 
pendently determined and it is assumed that selection effects are present. The model of 
figure 5 is remarkably close to the model found by a different method of model determina- 
tion presented in reference 10. The results obtained for a model for the universe when 
C is independently determined confirms the possible presence of selection in the obser- 
vational data and that the assumption that the parameter C is the same for all sources 
in the group being studied should not be used. 

NEW DATA 

Subsequent to completion of research for this report in mid 1968, new data on the 
age of the Galaxy have become available. Dicke (ref. 11) has published new ages based 
on nuclear cosmochronology, and new information relative to stellar ages, some not 
formally published, has been made available. A review of all information relative to the 
age of the Galaxy indicates that this age should now be considered to be greater than 
7.0 x 10^ years and possibly greater than 7.5 x 10^ years. These new lower boundaries 
for the age of the Galaxy give a much better fit with the time since the beginning of expan- 
sion for the model presented in figure 5 for which a D = 0.217 and q Q = 0.402. 

CONCLUDING REMARKS 

The use of red-shift — magnitude relation to analyze observed red shift and apparent 
magnitudes has, in the past, been predicated on the assumption that all sources have the 
same value of the parameter C. The models obtained by this method have values of the 
density parameter and the time since the beginning of expansion which were not consistent 
with the values of these parameters obtained by independent methods based on observation 
and models of the universe obtained by other means. Generally speaking, the density 
parameter was too high and the time since the beginning of expansion too short when com- 
pared with the independently determined values. 

When the effects of evolution and selection in the data were accounted for by an 
independent method for determining the parameter C (this method allows C to vary 
from source to source), more acceptable models for the universe were found. It was 
concluded that the assumption that all sources in the data have the same value of the 
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parameter C should not be used when analyzing observational data with the red-shift — 
magnitude relation. However, if the parameter C can be determined independently for 
each source, the red-shift — magnitude relation can be used for the analysis of observa- 
tional data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 10, 1969, 

188-41-01-01-23. 
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SYMBOLS 


Ag age of Galaxy 

C = 5 log 10 + M - 5 
c speed of light in a vacuum 

H 0 Hubble parameter 

i refers to members of group, i = 1, . . n 

K red-shift correction 

k curvature constant 

M absolute magnitude 

M rate of change of absolute magnitude 

AM incremental absolute magnitude due to evolution 

m apparent magnitude 

m c corrected observed apparent magnitude 

m cal calculated apparent magnitude 

m-Q calculated apparent magnitude for a Q = q Q = 1.0 

Am incremental apparent magnitude 

m mean apparent magnitude 

in^ mean apparent magnitude of kth group 

m^ mean corrected apparent magnitude of kth group 
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number of sources in group 
deceleration parameter 

scale factor that describes geometric history of universe 

present value of scale factor 

function that depends on curvature of space 

time since beginning of expansion 

time 

present time 
red shift 

c\ 

cosmical constant, 3 Hq (cr 0 - q 0 ) 
variance 


density parameter 
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APPENDIX B 

OBSERVATIONAL DATA USED IN THE STUDY 

The data used in the present paper with the red-shift — magnitude relation are the 
corrected apparent magnitudes and red shifts of galaxies that are members of clusters. 
Two principal groups of data fall into this category: data for 18 clusters reported by 
Humason, Mayall, and Sandage (ref. 8) and data for eight clusters observed by Baum. 

The red shifts and apparent magnitudes of the clusters of galaxies reported by Humason, 
Mayall, and Sandage (ref. 8) are as follows: 


Cluster 

Red shift, 
z 

Apparent 
magnitude , 
m 

Virgo 

0.004 

9.16 

0316 + 4121 

.018 

12.51 

1257 + 2812 

.022 

12.84 

1603 + 1755 

.036 

14.12 

2308 + 0720 

.043 

14.78 

2322 + 1425 

.044 

15.04 

1145 + 5559 

.052 

15.71 

0106 - 1536 

.053 

15.21 

1024 + 1039 

.065 

15.88 

1239 + 1852 

.072 

15.22 

1520 + 2754 

.072 

15.93 

0705 + 3506 

.078 

16.26 

1431 + 3146 

.131 

17.31 

1055 + 5702 

.134 

17.31 

0025 + 2223 

.159 

17.39 

0138 + 1840 

.173 

17.16 

0925 + 2044 

.192 

17.54 

0855 + 0321 

.202 

17.84 


The apparent magnitudes are photographic magnitudes and have been corrected for the 
aperture effect and the K correction by the methods given in reference 8. These inves- 
tigators give no uncertainties for these magnitudes; however, De Vaucouleurs (ref. 12) 
indicates that good photographic magnitudes carry an error from ±0.1 to ±0.2. It was 
assumed that these magnitudes carried a similar error. A survey of reports on red- 
shift measurements indicates that for a given observing system the uncertainty is in the 
fourth decimal place. However, when the red shift for a given source is determined by 
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more than one observing system, a comparison of the results indicates an uncertainty 
in the third decimal place. This uncertainty would seem to make the red shift of the 
Virgo cluster suspect since z = 0.004. However, the red-shift uncertainty mentioned 
applies to single spectrograms and the red shift of the Virgo cluster given is the mean 
for 73 members of the Virgo cluster; therefore, the red shift is probably accurate to 
three decimal places. 

Baum has never published a definitive paper on the magnitudes of the eight clusters 
he observed. In reference 13, Baum gives the red shifts he measured and the band 
passes of the photometer. Baum's data used in this study are those published by 
McVittie (ref. 2, p. 190). The red shifts and apparent magnitudes for the eight clusters 
of galaxies observed by Baum are as follows: 


Cluster 

Red shift, 
z 

Apparent 
magnitude , 
m 

Virgo 

0.004 

9.2 

1257 + 2812 

.022 

12.8 

1520 + 2754 

.072 

15.6 

1055 + 5702 

.134 

16.9 

0925 + 2044 

.192 

17.4 

0024 + 1654 

.290 

18.7 

1448 + 2617 

.350 

18.5 

1410 + 5224 

.460 

19.3 


Baum's magnitudes were obtained by use of an eight-color photoelectric photometer and 
no red-shift correction is needed. Again, no uncertainties are given but De Vaucouleurs 
(ref. 12) gives the uncertainty for good photoelectric magnitudes as ±0.02 to ±0.03. Inas- 
much as these uncertainties compare very well with those given for quasi-stellar objects, 
they were assumed to be applicable to Baum's data. The red-shift uncertainty discus- 
sion for the data of Humason, Mayall, and Sandage (ref. 8) applies to the data of Baum 
presented in reference 2. 

The smaller errors that appear to apply to Baum's magnitudes make it the better 
of the two sets of data; however, there are too few points to make an adequate analysis. 
The data of Humason, Mayall, and Sandage, although more extensive, cut off at a red shift 
that is too low for an adequate calculation of a model for the universe. 

These two sets of data when combined would constitute a more adequate set of data 
than either alone. The two sets of data were combined for this study as follows. The 
mean difference in apparent magnitude for the five sources common to each list was 
determined. Baum's data were then made fainter by the mean difference which was 0.18. 
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For common sources the average was taken; otherwise, the data of Humason, Mayall, 

and Sandage or modified Baum data were used. The combined data are presented in 

table I. Inasmuch as the method reduces Baum's data to those of reference 8, it was 

assumed that all data now represented corrected apparent photographic magnitudes and 

mm 

that the uncertainty was 0.10 to 0.20. 

In addition to the data for the red shift and apparent magnitude, an estimate of the 
age of the Galaxy and an estimate of the present density of the universe were also used. 
The age of the Galaxy can be determined from nuclear cosmochronology as was done by 
Fowler (ref. 14) and Hoyle (ref. 15) or from the age of the old population II stars in 
globular clusters. Fowler and Hoyle estimate the age of the Galaxy to be from 
1.5 + 0.5 x 10 10 years to 1.5 - 0.3 x 10-*-® years, an age that agrees well with the age of 
the population II stars in globular clusters. An age of the Galaxy of 1.5 x 10 10 years has 
been used in previous cosriiological studies (for example, ref. 16) and this value was 
adopted for this study. The present density of the universe is based on matter that 
shows luminous and dynamical properties. Abell (ref. 17) estimates the density on this 
basis and found that the density could range from 10 _ ^9 g-cm“3 to 10“31 g-cm - ^ with 
the most probable value being 10”30 g-cm - ^. The corresponding o Q interval is 
0.0045 = (J Q = 0.45. The value of cr 0 corresponding to most probable density of Abell is 
0.045. The methods used to estimate density do not rigorously account for intergalactic 
matter. There are two sets of observations in which intergalactic matter may have been 
detected. One observation by Koehler and Robinson (ref. 18) of absorption in the 21-cm 
line of hydrogen was interpreted as indicating an intergalactic medium with a density 
approximately 7 X 10"31 g-crn"^. The second observation by Henry and coinvestigators 
(ref. 19) of radiation in the 40 A to 60 A band was interpreted as indicating a hot interga- 
lactic plasma with a density of about 10 _ 5 proton-cm - ^ jf Koehler and Robinson are 
correct concerning intergalactic density, the value of op would not appreciably change. 
However, if Henry and his coworkers are correct about the intergalactic medium being a 
dense (10“5 proton-cm - 3) hot plasma, then o Q should be multiplied by 16. However, 
Henry and his coworkers have based these results on very limited observational data and 
until these results are confirmed, a value of o Q in the interval 0.45 S o 0 i 0.0045 is 
acceptable. This interval is centered on Abell's most probable value of the density. 

The idea that the universe was once in a highly condensed state which took the form 
of a primeval fireball has been espoused by many investigators. (See ref. 20 for a brief 
review of the idea.) One of the more interesting results to come from the primeval fire- 
ball idea is that a relic of the fireball survives as universal thermal radiation of about 
5° K to 6° K (ref. 21). A thermal radiation of about 3° K was detected by Penzias and 
Wilson (ref. 22) and has since been found to be isotropically distributed in the universe. 
The observational confirmation of the universal thermal radiation means that any model 
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which is proposed for the observed universe must be such that this radiation can be 
explained by one or more features of the model. Such a feature would be an extremely 
small radius before the start of expansion. 


20 



APPENDIX C 


A POSSIBLE METHOD FOR DETERMINATION OF THE PARAMETER C 

The determination of the parameter C for sources in any group of data permits 
both evolution and selection to be taken into account. The following method, used to obtain 
values of C in table HI, depends on two factors: a good estimate of the upper and lower 
bounds for ct 0 and q Q and a knowledge of the uncertainty in the measurement of appar- 
ent magnitude. The first step is to select a source between z = 0.0133 and z = 0.02 
and determine a value of C from the equation 

Ci = m - 5 log 10 z (Cl) 

In this red-shift interval, a Q and q Q being different from 1 do not introduce significant 
differences in the value of C. Let CT o^ > °o u anc * be the lower and upper bounds 

for cr 0 and q Q . Use the value of C determined from equation (Cl) (that is, C x ) and 
a °l * n e O ua tion (7) with M = 0 to compute an apparent magnitude m^ for each 

source in the data group. Repeat the procedure for a o u > c lo u to obtain m u . Next, 
evaluate the equation 

m n = 5 log 1Q z + Cj (C2) 

for each red shift and compute 

m Z x = m ll " 6 ( C3 ) 

and 

m u 1 = m n + e 

where e is the uncertainty in the magnitude measurement. Then, determine the red 
shift at which m^ = m^ and m u ^ = m u . Call these red shifts z^ and z u . The rela- 
tionship of the quantities used to determine C is shown in figure 6. The following set 
of inequalities is used to estimate C: 


For 

and 

the parameter C is 

N N 

IIA IIA 

& 

(mjj + ne) s m g (m ui + ne) 

Ci + ne 

z =% 
Z > Zu 

(m^ + ne) S m S (ra u + ne) 

Ci + ne 

Z > Z Z 
Z g Zu 

i + ne) § m £ ^m u ^ + ne) 

Ci + ne 

z > z, 

(mj + ne) £ m £ (m u + ne) 



Ci + ne 

Z > Zu 
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In the foregoing inequalities, n is zero or a positive or negative even integer. Posi- 
tive n represents decreasing brightness and negative n represents increasing bright- 
ness. The process for determining the possible values of C is to start with n = 0 and 
decrease n until all values of m are greater than m u + ne . The process is repeated 
for increasing n until all values of m are less than m^ + ne . If during this process 
m is equal to one of the boundary values when z is less than one of the critical values, 
that galaxy has two values of C. For example, if n = 0 and m = m^ and because 
m^ = m u - 2e, the galaxy can have values of C equal to C\ and Ci - 2e. Because 
for z > z u and z > z ^ the difference (m u + ne) - + ne) increases rapidly and if e 

is small, the distant sources can have a great many possible values of C. This large 
number of possible values of C can lead to a huge computing program of model deter- 
mination because all possible combinations of C must be considered. Both selection 
and evolution operate to make the distant galaxies brighter than the nearer ones. There- 
fore, the computing time can be reduced if only those combinations are considered where 
C for sources that have a red shift greater than one or both of the critical values of red 
shift is less than or equal to C for sources where red shift is less than the critical 
values. 


l0 3,nZ 



-m u boundary a^.q^ 


Figure 6.- Relationship of quantities used in the determination of the parameter C. 


This method for determining possible values of the parameter C is easily pro- 
gramed for digital computers and was used to determine the values of C in table III. 
The list of values of the parameter C given in table III is only one of many possible 
combinations. Once C has been determined, it can be used directly in equation (7) with 
AM =0 or it can be used to generate Am by subtracting 5 log^Q z + C from the 
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observed magnitude for use in equation (10) as was done in this report. Because some 
combinations of C can produce a model with negative density, the model determination 
program should be instructed to report only those combinations of C that give positive 
values of ct 0 . 
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